Response to the comments by Bernaschi et al. on 'Optimisation and parallelisation strategies for Monte Carlo simulation of HIV infection'

In recent years, the study of immune response behaviour through mathematical and computational models has attracted considerable efforts. The dynamics of key cell types, and their interactions, has been a primary focus in terms of building a picture of how the immune system responds to a threat. Discrete methods, based on lattice Monte-Carlo (MC) models, with their flexibility and relative simplicity have previously been used to model the immune system behaviour. However, due to speed and memory constraints, large-scale simulations cannot be done on a single computer. Key issues in the reduction of simulation time are code optimisation and code parallelisation. In this paper, optimisation and parallelisation solutions are discussed, with reference to existing MC simulation code for dynamics of HIV infection.